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This  is  one  of  three  reports  describing  the  Emissions  and  Dispersion  Modeling 
System  (EDMS).  All  reports  use  the  same  main  title — *  MICROCOMPUTER  MuuEL  FOR 
CIVILIAN  AIRPORTS  AND  AIR  FORCE  BASES— but  different  subtitles.  The  subtitles 
are : 


(1)  USER'S  GUIDE  -  ISSUE  2  -  (FAA-EE-88-3/ESL-TR-88-54) 

(2)  MODEL  DESCRIPTION  -  (FAA-EE-88-4/ESL-TR-88-53) 

(3)  MODEL  APPLICATION  AND  BACKGROUND  -  (FAA-EE-88-5/ESL-TR-88-55) 

The  first  and  second  reports  above  describe  the  EDMS  model  and  provide 
instructions  for  its  use.  This  is  the  third  report.  It  consists  of  an 
accumulation  of  five  key  documents  describing  the  development  and  use  of  the 
EDMS  model.  One  of  the  documents  shows  the  application  of  EDMS  to  the 
assessment  of  air  pollution  at  Stapleton  International  Airport. 

This  report  is  prepared  in  accordance  with  discussions  with  the  EPA  and 
requirements  outlined  in  the  March  27,  1980  Federal  Register  for  submitting  air 
quality  models  to  the  EPA. 
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INTRODUCTION 

Key  reports,  summarizing  the  Emissions  and  Dispersion  Modeling  System  (EDMS), 
have  been  assembled  in  order  to  describe  the  development  and  use  of  the  EDMS 
model.  This  document  is  the  repository  for  these  key  reports. 

These  reports  were  prepared  in  accordance  with  discussions  with  the 
Environmental  Protection  Agency  (EPA)  and  instructions  provided  in  the  Federal 
Register  (March  27,  1980)  for  the  submission  of  air  quality  models  to  the  EPA 
for  consideration  as  Guideline  models. 

MODEL  DEVELOPMENT  AND  USE  -  ITS  CHRONOLOGY  AND  REPORTS 

The  need  for  an  effective  model  to  evaluate  pollution  from  aircraft  was  first 
identified  during  the  analysis  of  Concorde  emissions  at  Dulles  International 
Airport  (IAD)  in  1976.  During  this  analysis,  it  was  noticed  that  existing 
models  can  not  accommodate  detailed  changes  in  power  setting,  speed,  and  ground 
track  as  aircraft  enter  different  operational  modes  at  an  airport. 

The  summary  report  of  the  analysis  of  Concorde  (and  other  aircraft)  emissions 
at  IAD  is  located  in  Appendix  A.  It  should  provide  insight  into  model 
development  problems  caused  by  the  unique  operational  characteristics  'f 
aircraft . 

The  need  for  a  flexible,  easy-to-use  model  to  analyze  aircraft  emissions  became 
more  apparent  in  the  late  1970's  when  both  the  Environmental  Protection  Agency 
(EPA)  and  the  International  Civil  Aviation  Organization  (ICAO)  started  to 
develop  the  rationale  for  the  imminent  engine  emission  standards.  The  basic 
question  to  be  answered  was,  "How  significant  a  source  of  pollution  are 
aircraf  t  ?" 

To  help  answer,  this  question  the  Federal  Aviation  Administration  (FAA) ,  the 
United  States  Air  Force  (USAF) ,  and  the  EPA  established  an  air  quality  study  to 
quantify  the  impact  of  aircraft  emissions  on  air  quality  at  airfields.  During 
this  study,  a  simple  model  was  developed  to  help  evaluate  monitoring  data. 
However,  this  model  was  quite  cumbersome  to  use  since  all  calculations  had  to 
be  made  by  hand.  A  summary  of  this  air  quality  study,  which  was  completed  in 
1980,  is  included  in  Appendix  B. 

The  nonavailabl ility  of  small  vet  powerful  computers  has  impeded  the 
development  of  a  simple,  flexible  model.  However,  in  the  late  1970?-, 
computers  having  this  capability  were  starting  to  become  available.  One  such 
computer,  the  Hewlett  Packard  97  (HP-97),  was  then  used  to  model  emissions  from 
an  aircraft  that  would  be  accelerating  down  a  runway  during  takeoff.  The 
resulting  model  is  described  in  Appendix  C. 

With  the  introduction  of  personal  computers  in  the  early  198u's,  the  HP-97  code 
was  reprogrammed  for  an  Apple  11+  microcomputer.  This  approach  led  to  the 
original  GIMM — Graphical  Input  Microcomputer  Model  which  was  completed  in 
1982.  Because  it  employed  the  more  powerful  Apple  computer,  GIMM  could  be 
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Realizing  the  effectiveness  of  GIMM  in  meeting  both  FAA  and  USAF  needs,  the  FAA 
and  the  USAF  issued  a  memorandum  of  understanding  (MOU)  to  formally  blend  the 
efforts  of  both  agencies.  This  MOU  documented  the  need  for  a  single  FAA/USAF 
microcomputer  model  to  evaluate  air  quality  at  both  airports  and  airbases. 

This  model — the  Emissions  and  Dispersion  Modeling  System  (EDMS) — incorporates 
the  emissions  and  dispersion  algorithms  of  the  original  GIMM  that  have  been 
speeded  up  and  processed  through  a  commercial  data  base.  EDMS  was  completed  in 
1985,  and  its  code  and  User's  Guide  were  released  to  the  general  public  in 
December  1985  as  report  FAA-EE-85-4 /ESL-TR-85-41 . 

Since  that  time,  major  modifications  have  been  made  to  the  original  EDMS  to 
enhance  its  usability  and  incorporate  an  integral  dispersion  model  into  its 
code.  A  prototype  of  this  expanded  model  was  completed  in  1986  and  was  used 
to  analyze  air  quality  at  Stapleton  International  Airport  in  conjunction  with 
the  building  of  a  new  runway  for  that  airport.  The  summary  report  for  this 
application  of  EDMS  is  contained  in  Appendix  E. 

Since  1986,  the  prototype  EDMS  has  been  incorporated  into  the  main  EDMS 
system,  and  the  final  model  has  been  submitted  to  the  EPA  as  an  agenda  item 
for  the  Fourth  Conference  on  Air  Quality  Modeling. 
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On  February  4.  1976.  the  Secretary  of  Transportation  ordered  the  FAA  to 
monitor  emissions  (and  noisei  at  Dulles  International  Airport  To  comply 
with  this  order,  it  was  necessary  to  measure  the  ambient  pollution  levels 
(background!  in  and  around  Dulles  Airport  and  to  trace  the  dispersion  of 
emissions  from  a  single  Concorde  aircraft.  While  the  more  conventional 
background  measurements  could  be  easily  performed,  there  was  no  known 
case  where  the  vertical  and  horizontal  profile  of  the  emission  plume  from 
a  single  aircraft  had  been  identified.  A  mobile  monitoring  program  was, 
therefore,  initiated  to  determine  if  the  emission  plume  of  a  taxiing  or  lak 
ing  off  aircraft  could  be  detected.  Special  instruments  were  required  to 
measure  the  discrete,  non-steady  state  nature  of  the  dispersion  of  the  air¬ 
craft  plume.  The  final  measurement  system,  which  consisted  of  contin¬ 
uously  recording  instruments  coupled  with  high-speed  chart  recorders, 
successfully  detected  emissions  from  a  single  aircraft.  Long  term  measure- 


On  February  4.  1976.  the  Secretary  of 
Transportation  ordered  the  FAA  to 
measure  the  noise  and  low-altitude 
emissions  of  the  Concorde  aircraft  in 
connection  with  a  16  month  trial  period 
for  that  aircraft  '  This  paper  describe-- 
the  low-altitude  emissions  portion  of  the 
program.  which  is  being  performed  ai 
Dulles  International  Airport 

Measurements  began  on  Mac  2V 
1976.  the  date  of  the  first  Concorde  de 
parture  from  Dulles  airport  Results 
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from  these  measurements  are  described 
through  Feb.  1977  ithe  cutoff  date  for 
preparation  of  this  paper  i  The  program 
schedule  is  listed  in  Tahle  I 

In  order  to  comph  with  the  Secre 
tary's  order  to  monitor  Concorde  emis 
sions.  it  was  necessary  to  measure  both 
the  change  in  ambient  air  quality  caused 
by  the  operation  of  a  single  Concorde 
aircraft  (single  event  t.  and  the  pollution 


background  on  and  off  the  airport 
While  the  more  conventional  back¬ 
ground  measurements  could  be  easily 
performed,  there  was  no  known  case 
where  the  vertical  and  horizontal  profile 
of  a  single  aircraft's  emission  plume  had 
been  identified  Past  studies-  A  had  been 
unsuccessful  in  such  identification,  be¬ 
cause  they  were  carried  out  in  high- 
background  congested  airports  where 


circuitous  taxiing  in  and  around  com¬ 
plex  terminal  areas  was  required. 

Background  levels  are  low  at  Dulles 
Airport,  where  most  aircraft  use  only 
one  isolated  ramp.  Aircraft  generally  do 
not  approach  the  terminal  where  pol¬ 
lution  from  numerous  sources  may  in¬ 
termingle.  A  mobile  monitoring  program 
was,  therefore,  started  to  determine  if 
the  emission  plume  from  a  taxiing  air¬ 
craft  could  be  detected.  Results  of  this 
program  showed  that  continuously  re¬ 
cording  instruments  coupled  with 
high-speed  recorders  could  detect 
emissions  from  a  single  aircraft.  Long 
term  measurements  of  background  and 
single  event  pollution  were  then 
begun 

Objective 

The  principal  objective  of  this  pro¬ 
gram  was  to  measure  the  effect  of  Con¬ 
corde  emissions  on  populated  areas  at  or 
near  Dulles  airport.  Air  quality  was  de 
termined  at  two  main  populated  areas, 
namely,  the  airport  itself  and  the  Ster¬ 
ling  Park  Community  These  locations 
are  shown  in  Figure  1 .  The  impact  of  the 
airport  (and  Concorde  I  emissions  on  the 
air  quality  at  Sterling  Park  was  deter¬ 
mined  by  measuring  the  pollution 
background  upwind  and  downwind  of 
the  airport.  The  impact  of  Concorde 
emissions  on  the  airport  itself  was  de 
fermined  by  measuring  the  change  in 
pollutant  concentration  levels  caused  by 
emissions  from  a  single  aircraft  as  it 
started,  taxied,  and  took  off.  The  dis¬ 
tance  from  the  taxiing  aircraft  source  at 
which  these  emissions  blend  into  the 
background  determined  the  effect  of 
Concorde  emissions  on  the  terminal 
area  Single  event  measurements  were 
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Figure  4.  Character ist»c  plot  of  NO,  coocentratKxi  diring 
takeotl 


minutes  Vertical  pollution  measure- 
ments  were  made  at  five  elevations  on 
two  vertical  towers. 

Site  Selection  and  Instrumentation 


Figure  3  TakeoM  station  locations 

also  made  to  determine  engine  emission 
rates  for  comparison  with  those  listed  in 
the  Concorde  environmental  impact 
statement  A  model  is  being  developed 
from  these  measurements  The  Con 
corde  influence'  area  show  n  in  Figure  1. 
will  be  determined  from  the  model 


for  taxi  and  start  idle  emissions  NO, 
was  the  tracer  for  takeoff  emissions 
Continuously  recording  instruments 
coupled  with  high-speed  chart  recorders 
were  employed  to  measure  plume  pas 
sage,  which  usually  lasted  less  than  two 


Site  selection  considerations  were:  1 
probable  success  in  detecting  an  event. 
2.  freedom  from  spurious  emissions;  3. 
frequency  at  which  aircraft  passed  the 
monitoring  sites;  4.  available  power;  5. 
wind  direction;  6.  noninterference  with 
other  airport  c-  -  iiions. 

Six  instrumented  trailers  plus  mobile 
equipment  were  moved  at  different 


Approach 

Six  t  -  quality  stations  and  two  ver 
tical  towers  were  employed  for  both 
background  and  single  event  mea-ure 
menus  Background  measurements  were 
performed  at  two  main  stations;  one 
measuring  the  airport  background  and 
the  other  measuring  the  community 
background  All  major  pollutants  < car 
bon  monoxide  ((’()>.  hydrocarbons 
tHfT.  nitrogen  oxides  iN(H.  nitrogen 
dioxide  i  NO  I.  ozone  If)  l  and  part  idesi. 
as  well  as  wind  direction  and  speed  were 
measured  continuously  at  these  loca¬ 
tions  A  third  station  i n it  tally  set  up  to 
monitor  takeofl  emissions  was  found  to 
be  useful  for  background  measurements 
This  station  recorded  CO.  NO,  and 
wind  speed  and  direction 

Single  event  measurements  were  used 
to  define  the  Concorde  influence  area 
and  to  provide  the  detailed  d.t'.i  (or 
background  pollution  anal- sis  Mea 
surements  were  made  at  three  stations 
and  on  two  towers  The  measurement 
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DECS  1926  TO  FEB.29  1977 
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Figure  6.  Characteristic  plot  of  CO  concentration  during 
taxi 


times  to  the  32  sites  shown  in  Figure  1 
for  pollution  assessment  purposes 
Location,  instrumentation,  measure 
ment  purpose,  and  operational  dates  are 
tabulated  in  Table  II.  Equipment  at  the 
following  sites  monitored  the  opera 
tional  modes  and  weather  parameters 
listed 

Site  Surnbcr 

Taxi  (surface)  4.  5,  10,  11, 24.  25. 

26 

Engine  Start  Idle  15.16.1? 

Takeoff  (precursor!  2.  3 

Takeoff  12.  13.  14 

Taxi  (plume  rise  20T.  2 IT.  22T. 

single  towerl  23T 

Taxi  I  plume  rise  27T.20T.21T. 

double  towerl  22T.  23T.  28TT 

29TT.  30TT. 
31TT.32TT 


Background  1,6 

Approach  19 

Climb  out  18 

Meteorological  8.  9,  6.  10.  12.  24 

(wind  speed  and 
direction) 

Meteorological  7 

(vertical 
temperature) 

Instrument  selection  was  influenced 
by  the  unique  nature  of  the  aircraft 
pollution  source.  Most  non  aircraft 
sources  are  steady  state  in  nature  and 
change  little  over  long  periods  of  time 
These  sources  are  amenable  to  long 
sampling  time  instrumentation  The 
emission  plume  from  a  moving  aircraft, 
however,  is  a  non-steady  state  puff  and 
undergoes  a  wide  concentration  excur 
sion  as  the  emission  plume  passes  over 


the  downwind  monitoring  station.  This 
passage  which  usually  takes  less  than 
two  minutes  requires  continuously  re¬ 
cording  instruments  and  high-speed 
chart  recorders  to  record  the  rapid  pas¬ 
sage  of  this  event.  This  equipment  was, 
therefore,  used  to  record  the  short- 
duration  passage  of  the  Concorde 
emission  plume. 

Equipment  for  background  and  single 
event  monitoring  are  listed  below: 

Carbon  Monoxide 

Intertech  Co. — URAS2 — NDIR 
Energetic  Sciences  Inc.,  Ecolyzer 
2600E 

Nitric  Oxide/Nitrogen  Dioxide 
Thermo  Electron  Co.  14B  Analys¬ 
er 

Monitor  Labs  Inc.,  8500  Calibra 

tor 

Total  Suspended  Particulates 
BGI-IIA  Hi  Volume  Sampler 
BGI-HCII  Standard  Calibrator 

Total  Hydrocarbons 

Beckman  Instruments  Inc..— 
Model  400 

Non-Methane  Hydrocarbons 
Beckman  Instruments  Inc.— 
Model  6800 

Ozone 

McMillan  Electronics  Co.— lido 
Analyser. 

1020  Ozone  Generator 

Wind  Speed  &  Direction 

Climet  Instruments  Co. — 01 1  - 1 
Wind  Speed  Transmitter.  (112-10 
Wind  Direction  Transmitter. 
060-10  Transmitter 


cowccmoi  aircra* T 


MO*  CONCORDE  AIRCRAFT 


•  MAY  is  TO  AUGUST  7  1*76 

•  AUGUST  I  TO  SEPTEMBER  K  1*7* 
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Figure  8.  Tower  installation 


Temperature 

Climet  Instrument  ( 'o. --015-3 
Temperature  Sensor.  06p. 10 

Translator 

A  major  consideration  in  site  selection 
was  the  traffic  pattern  hemp  u-ed  at  the 
airport  Most  commercial  aircraft  do  not 
operate  in  the  vicinity  o(  the  terminal, 
but  rather,  position  themselves  at  the  jet 
ramp  which  is  located  2d 00  ft  south  ot 
the  terminal  Airplanes  move  around 
this  ramp  in  a  clockwise  direction  For 
south  wind  operations,  which  are  pre 
dominant  during  the  summer  months, 
tin  airplane'  usually  proceed  Irom  the 
ramp  to  takeofl  runwav  1;  left,  which 
assures  the  shortest  possible  taxi  dis 
tanee  For  north  wind  operations,  air 
planes  proceed  from  the  ramp  to  runway 
1-left  Considering  this  traffic  pattern, 
the  most  effective  location  for  taxi 
monitoring  during  the  summer  month' 
is  at  the  turf  area  just  oft  the  northea' 
tern  edge  of  the  taxi  ramp  I  Figure  2l. 
M  onitoring  started  at  two  locations  i4 
and  5  i  A  t  hird  local  ion  was  added  first 
at  site  in  and  then  at  site  11*  provide 
the  three  points  needed  to  determine 
emission  dispersion  As  the  wind  shitted 
to  the  north  in  the  winter  month',  the 
three  taxi  monitoring  station'  were 
moved  to  the  other  side  of  the  taxi  ramp 
Mea-urement'  to  record  the  vertical 
pollution  profile  were  performed  at  t ho 
location  Power  for  all  taxi  monitoring 
was  provided  hv  an  FAA  15KW  Diesel 
electric  generator 

Location  of  the  takeoff  monitoring 
site'  <12.  13  and  Ml  was  determined 
through  analysis  of  precursor  measure 
ments  taken  at  sites  '2  and  It  Site-  12,  13 


spot-checked  at  Sites  15,  16.  and  17. 
Monitoring  for  queuing  was  initially 
planned  but  then  dropped  because 
queuing  did  not  occur  at  the  time  of 
Concorde  departure. 

Results 

Between  May  1976.  and  Feb  1977. 
the  Concorde  monitoring  system  re¬ 
corded  the  pollution  background  on  and 
off  the  airport,  and  emissions  from  air¬ 
craft  single  events  during  engine  start 
idle.  taxi,  and  takeoff.  Major  emphasis 
has  been  on  monitoring  the  pollution 
background  and  the  emission  plume 
transport  from  a  taxiing  or  taking  off 
aircraft 

Background  Measurements 

All  major  pollutants  were  measured 
at  Site'  1  and  6  The  background  was 
measured  to  relate  air  quality  on  and  off 
the  airport  and  to  compare  pollution  in 
the  vicinity  of  the  airport  with  the  na 
tional  ambient  air  quality  standards. 
These  data  are  reported  in  References 
6  9  and  the  analysis  will  be  described  in 
the  final  Concorde  monitoring  report. 

Takaoff  Measurements 

Emissions  at  Sites  12.  13,  and  14  were 
measured  from  Nov  1976toF'eb  1977. 
Site  locations  are  shown  in  Figure  3  and 
measurements  were  taken  al  the  three 
downwind  locations  A  characteristic 
trace  of  the  pollution  time  history  is 
shown  in  Figure  4.  Cumulative  peak 
concentrations  at  different  distances 
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Taxi  Measurements  (Surface) 

Emissions  at  sites  4,  5  and  11  were 
measured  from  May  to  Sept.  1976.  A 
characteristic  plot  of  air  quality  during 
Concorde  plume  passage  is  shown  in 
Figure  6.  Cumulative  plots  of  Concorde 
and  non-Concorde  emissions  are  shown 
in  Figure  7.  These  data  support  the  fol¬ 
lowing  trends. 

1.  The  average  peak  CO  concentration 
for  Concorde  is  1.7  times  higher  than 
the  average  concentration  of  the 
other  aircraft  monitored  at  a  location 
200  ft  downwind  from  the  taxiing 
aircraft 

2.  Emissions  from  Concorde  (and  other 
airplanes!  disperse  to  background 
levels  before  they  reach  the  terminal 
(2300  ft  from  the  ramp  taxiwayl 

3.  The  contribution  of  one  taxiing 
Concorde  to  the  hourly  average  CO 
concentration  of  all  other  sources  is 
less  than  0.1  parts  per  million  (ppm) 
at  locations  as  close  as  200  ft  from  the 
center  line  of  the  taxiway. 

Taxi  Meaauramants — (Tower) 

Two  tower  tests  were  performed  a 
single-tower  test  to  determine  vertical 
plume  characteristics  and  a  double- 
tower  test  to  determine  change  in  plume 
characteristics  between  the  two  tow 
ers 

The  single  tower  test  wee  performed 
on  Nov  ]  through  the  15th.  using  a  58  ft 
tower  with  lour  vertical  pollution  intake 
positions  The  second  test  was  started 
on  Feb  20.  1977,  for  a  planned  week 
time  period  and  consisted  of  the  first  58 
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emissions  will  be  reported  after  the  1 
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reduce  to  levels  undetectable  from  the 
background  within  2000  ft  of  the  taxiing 
aircraft.  Concorde  (and  other  aircraft) 
single-event  emissions  contribute  less 
than  0.1  ppm  of  CO  to  ambient  air  con¬ 
centrations  at  locations  as  close  as  200  ft 
from  a  taxiing  aircraft,  when  averaged 
over  a  one  hour  time  period. 

Tower  measurements  show  that  the 
hot  emission  plume  tends  to  lie  close  to 
the  ground  and  does  not  rise  signifi¬ 
cantly  at  monitoring  station  locations. 
Specific  relationships  between  surface 
and  higher-level  concentrations  will  be 
evaluated  later  on  in  the  program. 


A  simple  model  is  being  developed 
from  the  measurement  data  The  model 
describes  the  downwind  concentrations 
of  the  emissions  from  a  single  aircraft 
passage  and  sums  up  the  emission  con 
trbutions  of  a  number  of  single  aircraft 
passage*-  over  longer  time  periods 
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APPENDIX  B 

THE  INFLUENCE  OF  AIRCRAFT  OPERATIONS  ON  AIR  QUALITY  AIRPORTS 


the  influence  of  aircraft  operations  on 
air  quality  at  airports 


Howard  Segal 

Federal  Aviation  Administration 
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Emission  standards  for  aircraft  engines  were  promul¬ 
gated  in  1973  after  it  was  determined  that  these  engines 
were  significant  sources  of  pollution  around  airports. 
Since  that  time,  new  information  has  become  available 
on  the  modeling  and  monitoring  of  aircraft  emissions 
and  in  March  1978  the  Environmental  Protection 
Agency  (EPA)  announced  its  intention  to  amend  the 
1973  standards.  Included  in  this  announcement  was  the 
establishment  of  a  joint  FAA/EPA  air  quality  study 
which  included  the  reassessment  of  the  impact  of  air¬ 
craft  emissions  on  air  quality  around  airports.  This 
paper  presents  the  results  of  this  study  which  includes 
the  assessment  of  air  quality  at  five  commercial  and  One 
general  aviation  airport.  Also  presented,  are  the  pre¬ 
liminary  results  of  research  performed  since  the  com¬ 
pletion  of  the  study.  It  is  concluded  that  aircraft  emis¬ 
sions  have  a  smaller  impact  on  air  quality  than  had  been 
estimated  in  studies  that  were  performed  prior  to  the 
promulgation  of  the  aircraft  engine  emission  standards 
in  1973. 


The  1970  amendments  to  the  Clean  Air  Act1  directed  the 
Environmental  Protection  Agency  (EPA)  to  establish  emis¬ 
sion  standards  for  aircraft  and  aircraft  engines,  if  such  emis¬ 
sions  are  judged  to  cause  or  are  likely  to  cause  or  contribute 
to  air  pollution  which  endangers  public  health  or  welfare.  The 
1970  amendments  also  directed  the  EPA  to  conduct  a  study 
of  the  extent  to  which  aircraft  emissions  affect  air  quality  in 
Air  Quality  Control  Regions  (AQCR)  throughout  the  United 
States.  Based  upon  information  available  in  the  early  1970s, 


Since  that  time,  major  advances  have  been  made  in  the 
techniques  for  monitoring  and  modeling  aircraft  emissions. 
On  March  24,  1978,  a  Notice  of  Proposed  Rulemaking 
(NPRM)  was  published  in  the  Federal  Register  to  announce 
the  intention  of  the  EPA  to  amend  the  1973  engine  emission 
standards.  Included  in  the  NPRM  was  the  establishment  of 
a  joint  Federal  Aviation  Administration  (FAA)ZEPA  air 
quality  study  to  relate  aircraft  emissions  to  ambient  air 
quality. 

In  setting  up  this  air  quality  study,  it  was  decided  first  to 
review  data  generated  both  before  and  after  the  promulgation 
of  the  engine  emission  standards  to  determine  the  complete¬ 
ness  of  these  data  in  establishing  an  air  quality  basis  for  the 
engine  emission  standards.  During  this  review,  it  was  found 
that  the  magnitude  of  the  initial  pollutant  dilution  caused  by 
exhaust  gas  heat  and  turbulence  was  only  first  measured  in 
1976.  Therefore,  modeling  results  prior  to  1976  could  be 
subject  to  substantial  error.  New  pollution  monitoring  pro¬ 
grams  were  then  initiated  to  get  additional  plume-related 
information.  The  resulting  monitoring  analysis  programs 
summarized  in  this  paper  involved  the  coordinated  efforts  of 
the  FAA,  the  EPA,  and  the  Air  Force,  Argonne  National 
Laboratory  (ANL)  and  Environmental  Research  and  Tech¬ 
nology,  Incorporated  (ERT).  A  detailed  technical  report  has 
also  been  issued.2 

Approach 

The  need  for  engine  emission  control  is  determined  through 
the  evaluation  of  economic,  technological,  and  air  quality  data. 
This  paper  addresses  the  evaluation  of  air  quality  data 
through  the  modeling  and  monitoring  of  aircraft  emissions  at 
a  number  of  airports.  Essential  to  this  evaluation  is  the  de¬ 
velopment  of  good  air  quality  data  for  model  validation.  But 
it  was  not  until  1976  that  the  emission  plume  from  aircraft  was 
successfully  isolated  from  other  airDort  sources.  This  nlunte 


of  the  Concorde  measurement  program,3  permitted  the  ver¬ 
tical  profile  and  trajectory  of  the  aircraft  emission  plume  to 
be  quantified  for  the  first  time.  A  typical  event  is  depicted  in 
Figure  1  and  an  ensemble  of  such  events,  spanning  a  wide 
variety  of  meteorological  conditions  and  aircraft  types  enabled 
development  of  a  plume  rise  equation  for  taxiing  aircraft  ex¬ 
haust  plumes. 

The  pollution  assessment  strategy,  which  is  described  in 
Figure  2,  was  directed  toward  quantifying  the  initial  size  and 
height  of  the  emission  plume  from  individual  aircraft  prior 
to  the  dispersion  of  this  plume  in  the  ambient  air.  This  in¬ 
formation  was  incorporated  into  the  Airport  Vicinity  Air 
Pollution  (AVAP)  model4  which  was  then  used  to  calculate 
pollutant  concentrations  at  the  three  large  airports  that  were 
evaluated  in  conjunction  with  the  1973  engine  emission 
standards. 

Monitoring  and  modeling  program  characteristics  at  four 
airports  are  listed  in  Table  I  Modeling  results  made  use  of  the 
following  "worst  case"  conditions  in  an  attempt  to  reflect  the 
implied  meteorological  and  activity  conditions  of  the 
NAAQS 

Averaging  time — I  hr 
Pasquill/Gifford  stability  class — E 
Wind  speed — 1  m/sec 
Aircraft  activity — Peak  levels 

Receptor  location— 750  m  downwind  from  the  runway/ 
taxiway.  (This  is  the  characteristic  distance  at  which  the 
general  public  might  first  be  exposed  to  pollution  from 
aircraft.) 

Pollutant  Considerations 

Four  pollutants  emitted  by  aircraft,  CO,  HC,  nitrogen  oxide 
(NO)  and  nitrogen  dioxide  (NOj)  (NO  +  NO2  is  referred  to 
as  NO, )  have  been  judged  in  the  past  to  be  significant.  Only 
NO2  and  CO  are  considered  in  detail  in  this  paper  because 
their  concentrations  can  be  directly  compared  to  an  appro¬ 
priate  NAAQS.  The  impact  of  aircraft  HC  and  NO  emissions 
on  oxidant  levels  was  not  addressed  in  this  study  because  of 
the  state-of-the-art  limitations  in  the  modeling  and  related 
monitoring  of  photochemical  reactions,  especially  when  the 
reactive  species  are  not  well  defined.  The  rationale  for  possibly 
controlling  these  precursor  pollutants  must  arise  from  other 
considerations. 


Flpar*  2.  Analysis  procedure 


Results 

Measurement  and  modeling  results  at  five  airports  are 
summarized  in  Figure  3.  Aircraft  related  concentrations  have 
been  determined  through:  statistical  analysis  of  measure¬ 
ments,  submodeling  of  aircraft  operations  at  airports,  and 
AVAP  modeling  of  aircraft  operations  at  airports.  While  there 
are  uncertainties  in  each  of  these  analysis  methods,  the  use 
of  three  independent  methods  permits  one  to  make  a  com¬ 
parison  of  the  consistency  of  the  three  results. 

The  right  hand  column  of  Figure  3  represents  the  results 
of  the  statistical  analysis  of  pollution  measurements  after  the 
background  had  been  subtracted.  The  values  selected  are 
based  either  on  data  extrapolation  to  reflect  the  one-hour  per 
year  that  the  short-term  standard  may  be  exceeded  or  else  on 
the  average  number  of  flights  per  hour  times  the  average  dose 
impact  per  flight  to  reflect  the  annual  standard.  The  middle 
column  represents  hourly  average  concentrations  that  were 
estimated  with  a  submodel  that  had  been  verified  with  met 
surement  data.  The  left  hand  column  represents  the  hourly 
average  concentrations  determined  with  the  AVAP  model, 
after  it  has  been  adjusted  to  reflect  measurements  of  initial 
exhaust  plume  size  and  height,  peak  aircraft  activity  levels, 
and  observed  aircraft  times  in  mode.  The  three  airports  used 
in  the  AVAP  assessment  were  the  same  airports  assessed  in 
conjunction  with  the  1973  engine  emission  standard:  John  F. 
Kennedy  (Jr  K),  O’Hare  International  (ORD),  and  Los  An¬ 
geles  International  (LAX). 

CO  Concentrations 

In  the  right  hand  column  of  Figure  3  it  is  seen  that  a  peak 
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aircraft  that  were  operating  at  Lakeland  Airport,  FL.  Sub¬ 
model  calculations  at  LAX  (middle  column)  and  AVAP  cal¬ 
culations  at  JFK,  ORD,  and  LAX  (left  column)  all  yield  peak 
hourly  CO  concentrations  in  the  range  of  4-7  ppm. 

From  all  these  data,  it  is  seen  that  no  estimate  of  peak 
hourly  average  CO  concentration,  at  a  source- receptor  dis¬ 
tance  of  750  m,  exceeds  one-fifth  of  the  NAAQS  for  CO. 

N03  Concentrations 

Measurement  and  model  results  have  been  compared  with 
standards,  taking  into  account  the  number  of  times  that  the 
standards  can  be  exceeded.  Two  types  of  standards  were 
considered:  a  long-term  NAAQS  (one  year  arithmetic  mean) 
and  a  possible  short-term  standard  (under  consideration). 
Since  one  would  expect  this  possible  short-term  standard  to 
reflect  the  characteristics  of  other  short-term  NAAQS,  this 
study  employs  concentration  averaging  times  of  1  hr  with  an 
expected  exceedance  of  one  time  per  year. 

NO2  data  have  been  analyzed  to  reflect  both  long  and 
short-term  standards.  With  regard  to  the  long-term  NAAQS, 
0.005  ppm  has  been  measured  at  300  m  from  the  source  when 
aircraft  depart  at  an  hourly  rate  of  10.1  airplanes/hr.  These 
data  are  shown  in  Figures  3  and  4.  Even  allowing  for  an  ad¬ 
ditional  factor  of  two  variation  resulting  from  seasonal  vari¬ 
ation  of  oxidant  levels  and/or  wind  direction,  such  levels,  while 
small  relative  to  the  long  term  NAAQS,  would  be  even  smaller 
at  the  750  m  source-receptor  distance  used  in  this  study. 

Study  data  were  also  processed  to  reflect  short  averaging 
times,  since  a  possible  short-term  standard'1  appeared  immi¬ 
nent  at  the  time  that  the  study  was  started.  As  of  the  writing 
of  this  paper,  the  possible  standard  had  not  materialized.  A 
different  criterion  was  therefore  sought.  In  1977,  a  World 
Health  Organization  (WHO)  task  group  selected  0.5  ppm  (1 
hr  average)  as  their  estimate  of  the  lowest  observed  health 


effect-level  for  short-term  exposure.*  This  value  is  plotted  in 
Figure  3.  (The  0.5  ppm  limiting  level  is  not  to  be  considered 
a  standard  which  would  have  to  include  an  adequate  margin 
of  safety.) 

The  decision  on  how  much  NO2  is  chargeable  to  aircraft  is 
difficult  to  make  since  NO2  is  generated  both  in  the  engine 
during  combustion  and  in  the  ambient  air  through  the  reaction 
of  ambient  pollutants,  including  ozone,  with  engine  produced 
NO  (approximately  95%  of  engine  NO,  is  released  as  NO). 
While  most  ambient  pollutants  react  quite  slowly  with  NO, 
03  reacts  very  quickly  with  NO  making  it  the  predominant 
precursor-NO  reaction  at  doee-in  locations  (less  than  1000  m 
between  source  and  receptor). 

Recent  measurements  of  the  O3-NO  reaction  in  aircraft 
plumes  at  O’Hare  International  Airport  (ORD)  have  sup¬ 
ported  the  “ozone  limiting”  approach  for  estimating  total  NO? 
concentrations  from  aircraft  at  cloae-in  locations.  Using  this 
approach,  the  NO2  values  listed  in  the  AV  AP  column  of  Figure 
3  were  determined  by  adding  the  NO2  produced  in  the  engine 
to  the  NO2  produced  by  the  reaction  of  engine  produced  NO 
with  ambient  O3.  Since  there  is  usually  a  surplus  of  NO  at 
critical  pollution  assessment  times,  03  concentration  is  the 
major  limit  to  the  amount  of  NO2  produced  in  the  ambient  air 
at  close-in  locations.  Since  high  63  levels  would  be  expected 
during  worst  case  conditions,  the  limiting  level  of  the  03 
NAAQS  (0.12  ppm)  was  assumed.  As  a  result,  0.12  ppm  of  Os 
would  react  with  an  equivalent  concentration  of  engine  pro¬ 
duced  NO  to  produce  0.12  ppm  of  N02. 

This  0.12  ppm,  when  added  to  the  0.02  ppm  produced  in  the 
engine  results  in  a  total  NO2  concentration  of  0.14  ppm.  This 
value  is  plotted  in  the  left  hand  column  of  Figure  3. 

The  NO2  value  of  0.2  ppm  in  the  right  column  represents 
the  average  of  the  NO2  values  determined  from  DCA  mea¬ 
surement  data.  Again,  as  with  the  CO  and  NO2  dsta  reported 
earlier,  concentrations  at  the  measurement  distance  of  300  m 
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shown  in  Figure  4  would  be  even  lower  at  the  characteristic 
750  m  source-receptor  distance  used  in  this  study. 

HC  Concentrations 

The  results  of  AVAP  modeling  at  JFK,  ORD,  and  LAX, 
tabulated  in  Figure  3,  indicate  that  pealt  hourly  concentrations 
from  aircraft  are  approximately  5  ppm.  Since  there  is  no  HC 
NAAQS,  this  concentration  cannot  be  related  to  any  partic¬ 
ular  standard. 

Conclusions 

Data  from  Figure  3  on  the  impact  of  aircraft  emissions  on 
air  quality  can  be  summarized  as  fillows: 

For  CO,  1  hr  average  concentrations  do  not  exceed  7  ppm 
when  concentrations  are  determined  under  conditions 
compatible  with  the  NAAQS.  This  value  is  small  relative 
to  the  35  ppm  limiting  level  of  the  NAAQS. 

For  NO2,  annua]  average  concentrations  are  only  10-20% 

Kli  inc  01 * _ Kin - .-*-■ —  -v  * 


The  impact  of  aircraft  HC  (and  NO, )  emissions  on  oxidant 
levels  are  not  addressed  in  this  study  because  of  the  state- 
of-the-art  limitations  in  the  photochemical  modeling  and 
monitoring  of  aircraft  emissions. 

CO  and  NO2  pollution  from  aircraft  appears  to  be  small 
relative  to  pollution  limits  in  the  NAAQS.  This  is  caused  by 
enhanced  initial  dispersion  due  to  the  heat  and  turbulence 


associated  with  j«t  exhaust  plumes  and  the  strong  localization 
of  aircraft  emissions  at  areas  near  the  ends  of  runways  and, 
consequently,  quite  remote  from  locations  of  public  exposure. 
These  factors,  not  present  in  the  case  of  automobile  pollution, 
act  to  mitigate  the  significance  of  pollution  from  aircraft. 
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INTRODUCTION 

Atmospheric  dispersion  models  are  mathematical  expressions  that  combine 
source  emissions  with  meteorological  parameters  to  produce  air  quality 
estimates  at  specified  receptor  locations.  At  airports  where  many 
sources  and  receptors  are  in'  olved,  refined  models  such  as  the  Airport 
Vicinity  Air  Pollution  model  (AVAP)(1)  are  used  to  determine  air  quality. 
However,  where  few  sources  and  receptors  are  involved,  screening  models 
are  very  attractive  for  identiying  the  need  for  further  analysis  with 
the  more  refined  models.  This  report  describes  one  of  these  screening 
models,  SIMPLEX  "A". 

This  report  describes  the  mathematical  basis  for  the  model,  lists  the 
progi im,  and  explains  the  steps  taken  to  compute  pollution  dosage. 

Special  program  features  are  described  and  two  sample  problems  are 
solved . 

The  experienced  user,  who  is  primarily  concerned  with  running  a  specific 
problem,  may  bypass  the  descriptive  sections  of  this  report  and  proceed 
directly  to  the  "Sample  Problem-Program  Operation"  section  on  page  5. 

MODEL  DESCRIPTION 


SIMPLEX  "A",  which  has  been  programmed  for  the  Hewlett  Packard  67  and  97 
desk  calculators,  addresses  emissions  during  takeoff.  Additional  SIMPLEX 
models  are  being  developed  to  determine  the  air  quality  impact  from 
taxiing  and  queueing  aircraft  as  well  as  from  ground  vehicles  at  the 
airport.  The  model  is  particularly  useful  at  small  airports  and  at 
those  airports  having  only  a  few  dominant  sources. 


SIMPLFX  "A"  uses  the  same  Gaussian  formulation  employed  in  many  of  the 
refined  models  listed  in  the  Environmental  Protection  Agency's  (EPA) 
guidelines  on  air  quality  models.  It  accomplishes  its  function  by 
simplifying  many  of  the  detailed  features  of  the  more  refined  models. 

The  model  is  an  integrated  puff  model  for  an  accelerating  point  source. 
Downwind  receptors  are  assumed  to  be  at  gro  "d  level  (z=0)  and  receive 
pollution  doses  from  each  emission  puff.  Figure  1  describes  the  source- 
receptor  geometry  where  the  dose  from  each  emission  puff  is  summed  at  a 
receptor  to  give  a  total  dose  due  to  a  complete  takeoff  event.  Concentrations 
are  measured  in  parts  per  million  (ppm)  of  Nitrogen  Dioxide  (NO^)  where 
the  complete  conversion  of  Nitrogen  Oxides  (NO  )  to  NO2  is  assumed.  In 
cases  where  Carbon  Monoxide  (CO)  concentrations  are  required,  NO2 
calculations  can  be  factored  appropriately.  The  total  dose  at  point 
x,y,0  is  given  by  the  equation: 
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SYMBOL  DEFINITION 


UNITS 


^ *  receptor  exposure  of  dose  ppm-sec.  (NO^) 

X  ■  downwind  distance  in  the  direction 

of  the  mean  wind  meters  (m) 

Y -  crosswind  dis'ance  m 
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H  = 


height  above  ground  level 
standard  deviation  of  plume  concen¬ 
tration  in  the  vertical  direction 
standard  deviation  of  plume  concen¬ 
tration  in  the  crosswind  direction 
wind  speed 

total  emissions  during  an  emission 
release 

effective  height  of  emissions 


m 

m 

m 

m/sec 

grams 

m 


The  program  is  printed  out  in  Figure  2.  The  registers,  labels,  inputs 
and  outputs  are  listed  in  Figures  3  and  4. 


SPECIAL  PROGRAM  FEATURES 


Standard  Deviations  of  Plume  Concentration  (sigma  (ff~) ) 

A  subprogram  was  employed  to  determine  the  standard  deviation  of  plume 
concentration  in  the  horizontal  (crosswind)  and  vertical  directions. 
This  subprogram  was  based  upon  the  assumption  that  pollution  disperses 
according  to  the  power  law  expression. 
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<T  =  Kxb 


or,  in  straight  line  form 
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Log  <J*=  b  Log  x  +  Log  K 


The  exponent  'b"  governs  the  rate  of  pollutant  dispersion  and  the  coefficient 
"K"  depends  upon  atmospheric  stability. 


Analysis  of  the  Pa squ i 11 /G i f f ord  curves*  (3)  used  in  most  dispersion 
models  shows  that  for  stability  classes  "B"  through  "E",  single  straight 
lines  are  approximated  when  (Tv  and CTz  values  are  plotted  logarithmically 
against  downwind  distance  up  to  a  source-receptor  distance  of  1000 
meters.  It  is  also  seen  that  these  straight  lines  have  the  same  slope. 


With  the  realization  that  C  as  a  function  of  four  stability  classes  can 
be  represented  as  single  straight  lines  with  the  same  slope  (0.9), 
equation  2  can  be  rewritten  as: 
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C Tv  =  KjX 

0~  7  =  k2x 


0.9 

0.9 


when  x  does  not  exceed  1000  meters 


The  values  of  K  and  which  are  listed  in  the  program  printout  (Figure 
2)  were  obtained  by  solving  equations  4  and  3  for  and  after 
substituting  values  for  x,^"y  an d(T*z. 

The  program  not  onlv  allows  for  calculation  of  the(Ty  and  tf"z  values  in 
equations  4  and  3  but  also  has  provisions  to  input  values  for  initial 
sigmas  (C o  and(T'v)  in  order  to  account  for  the  enhanced  dispersion 
- — — Aral  nrUv  ±pf  pxhfluqt  .  (Thi«;  enhanrpd  disnersion 
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Because  of  technical  difficulties,  it  has  not  been  possible  to  determine 

CTo  values  from  measurements  taken  at  airports  during  high  thrust  airplane 

takeoff.  However,  plume  measurements  have  been  made  during  low  thrust 

operations  at  Dulles  (4)  and  Los  Angeles  International  (7)  airports. 

Average  values  from  measurements  taken  at  these  two  airports  (8  meters 

for  (To  and  16  meters  for C7o  )  are  incorporated  in  the  model. 
z  y 

Plume  Height 

Because  of  the  lack  of  experimental  data  to  support  plume  rise  theories 
for  taking-off  aircraft,  special  plume  rise  algorithms  have  not  been 
incorporated  into  the  model.  While  research  is  planned  in  this  area, 
until  this  research  is  completed,  it  was  assumed  that  the  plume  height 
was  at  least  as  high  as  the  airplane  engines.  An  average  value  for  this 
parameter  is  four  meters  for  airplanes  operatng  at  a  typical  large 
airport. 

Stability  Class 

Pasquill/Gif ford  stability  classes  "B",  "C",  "D",  and  "E"  are  expected 
to  prevail  at  the  airport  during  the  times  of  air  quality  assessment. 
Turner  (3)  gives  a  detailed  description  of  the  characteristics  of  each 
stability  class.  A  particular  stability  class  is  identified  bv  a  range 
of  wind  speed,  solar  radiation  intensity,  and  cloud  cover.  Values  for 
these  parameters  can  be  obtained  from  local  National  Weather  Service  or 
observer  personnel. 

Winds 

The  coordinate  system  is  oriented  to  the  runway  on  which  the  aircraft 
are  assumed  to  be  operating.  Since  aircraft  usually  take  off  into  the 
wind,  wind  angles  are  measured  only  from  0  to  90  degrees  on  either  side 
of  the  runway.  For  example,  a  zero  degree  wind  would  blow  directly  down 
the  runway;  a  90  degree  wind  would  blow  perpendicular  to  the  runway. 

Acceleration 

Aircraft  performance  manuals  can  be  used  to  determine  acceleration 
during  takeoff.  However,  the  program  has  been  structured  to  accept  an 
average  takeoff  acceleration  and  performance  information  should  be 
adjusted  to  average  acceleration  values. 

Emission  Tail 


During  the  operation  of  a  jet  engine,  the  high  velocity  of  its  exhaust 
gases  creates  an  emission  tail  which  can  extend  for  a  considerable 
distance  behind  the  aircraft  (Figure  5).  This  tail  is  simulated  by- 
assuming  a  value  for  its  length  and  by  assuming  a  finite  number  of 
points  along  the  tail  at  which  emissions  are  considered  to  be  released. 


i  f  ► 


»  rt  A  ►  V\  j* 
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The  program  assumes  a  225  meter  emission  tail  with  three  emission  release 
points  located  75  meters  apart  in  the  tail  (see  Figure  5).  The  model  is 
programmed  to  index  the  emission  starting  point  75  meters  further  down 
the  tail  after  each  iteration  sequence  is  completed.  The  first  and  last 
points  in  the  tail  are  37.5  meters  from  the  ends  of  the  tail.  (The  tail 
starts  at  the  exit  plane  of  the  engine.) 

Vertical  Dispersion  Lid 


Calculations  under  a  variety  of  assessment  conditions  showed  that  a  lid 
on  vertical  dispersion,  i.e.  an  inversion  "cap",  had  an  insignificant 
effect  on  concentration  at  the  short  downwind  distances  (less  than  1000 
meters)  employed  in  assessing  aircraft  pollutant  impact.  An  algorithm 
to  account  for  this  phenominon  is,  therefore,  not  included  in  the  program. 

Iteration  Interval 

From  past  program  use,  a  one-second  iteration  time  is  recommended. 

Using  this  iteration  time  interval,  the  dose  calculation  can  be  completed 
in  less  than  15  minutes  at  a  source  receptor  distance  of  300  meters. 

Dosage  Output 


Total  dosage  is  printed  out  at  the  end  of  each  iteration  sequence.  The 
program  is  stopped  when  the  dosage  reaches  a  maximum  value.  Output 
units  are  parts  per  million-seconds  (ppm-seconds) .  To  determine  the 
average  concentration  over  a  one-hour  time  period  (for  comparability 
with  a  particular  short  term  standard)  the  dosage  must  be  divided  by 
3600  seconds. 

SAMPLE  PROBLEM 

The  step  by  step  procedure  for  solving  the  sample  problems  is  described 
in  this  section.  While  this  procedure  is  structured  for  a  single  aircraft, 
the  same  procedure  can  be  used  for  any  number  of  aircraft  by  treating 
them  as  one  large  aircraft. 

Preparation  of  Data  For  Program  Execution 

From  Figure  1  it  is  seen  that  the  Case  1  receptor  is  located  337.5 
meters  downwind  from  the  aircraft  (as  measured  along  the  runway)  and  200 
meters  abeam  of  the  runway  centerline.  The  Case  2  receptor  is  located 
262,5  meters  upwind  of  the  aircraft  and  100  meters  abeam  to  it. 

The  objective  of  this  problem  is  to  determine  the  air  quality  impact  of 
747  N0X  emissions  (reported  as  NO2)  during  takeoff.  During  this  takeoff, 
it  was  assumed  that  a  5-meter  per  second  wind  was  blowing  at  30  degrees 
to  the  runway  centerline  and  that  Pasquill/Gif ford  stability  class  "E" 
prevailed.  The  747  was  assumed  to  have  a  constant  takeoff  acceleration 
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The  following  procedure  was  used  in  solving  the  problems: 

Source  emissions  were  obtained  from  AP-42  supplement  10  (8),  where  747 
NO  emissions  are  listed  at  215.3  kilograms  per  hour  per  engine  or  60 
grams  per  second  per  engine.  Since  the  747  has  four  engines,  the  total 
emission  rate  was  240  grams  per  second.  To  accommodate  the  three  emission 
release  points  in  the  "tail"  (see  Figure  5)  this  rate  was  divided  by  3 
to  reduce  its  value  to  80  grams  per  second  per  "tail"  release  point. 
Selecting  an  iteration  time  of  one  second  and  multiplying  it  by  the 
emission  rate  results  in  the  release  of  80  grams  of  N0^  per  puff. 

Xfbz  and  <Toy  values  of  8  and  16  meters  respectively  were  selected  from 
the  Standard  Deviation  of  Plume  Concentration  section  of  this  report  and 
a  plume  height  of  4  meters  was  selected  from  the  plume  height  section. 

The  beginning  time  was  set  at  zero  by  inputting  the  iteration  time  (one 
second)  and  assigning  a  negative  sign  to  it.  The  Case  1  receptor  is 
downwind  of  the  aircraft  giving  it  a  positive  sign  (see  Figure  1).  The 
Case  2  receptor  is  upwind  of  the  aircraft  giving  it  a  negative  sign. 

The  airplane  to  receptor  distance  is  converted  to  a  "tail"  to  Teceptor 
distance  (at  the  firs'"  "tail"  emission  point,  see  Figure  5)  by  subtracting 
37.5  meters  from  the  former  to  uniformly  space  the  three  tail  release 
points  over  the  225  meter  tail  length.  The  resulting  distance  between 
the  receptors  and  the  first  point  in  the  emission  tail  is  +300  meters 
for  Case  1  and  -300  meters  for  Case  2. 

Program  Operation 

Load  the  Program 

Before  loading  the  program  the  "on-off"  switch  should  be  in  the  "on" 
position  and  the  "run-program"  switch  should  be  in  the  "run"  position 
(for  the  HP-97  the  "trace-manual-norm"  switch  should  be  in  the  "manual" 
position).  The  program  can  then  be  loaded  into  the  calculator  by  first 
pushing  the  number  1  end  of  the  magnetic  tape  strip*  into  the  slot  in 
the  upper  left  hand  portion  of  the  HP-97  calculator.  (On  the  HP-67 
calculator  the  slot  is  located  on  the  right  hand  side.)  When  the  strip 
comes  out  the  other  side,  turn  it  around  to  the  number  2  end  and  push  it 
through  the  slot  a  second  time.  The  program  is  now  loaded  into  the 
calculator  and  the  tape  strip  which  has  come  out  the  back  of  the  calculator 
can  be  stored  in  the  horizontal  slot  just  under  the  calculator  switches. 

Input  Data 

Inputs  for  the  Case  1  and  Case  2  problems  are  listed  in  Figures  6  and  7 
and  a  printout  of  the  results  is  listed  in  Figures  8  and  9.  The  Case  1 
problem  is  solved  by  first  entering  the  values  for  the  six  input  parameters 
listed  in  Figure  6  into  the  Primary  Register  by  depressing  the  following 
keys:  80  STO  0  8  STO  2  4  STO  B  30  STO  D  1  STO  E  5  STO  I.  Any  input 
errors  can  be  erased  by  depressing  the  CLx  key  or  by  turning  off  the 
calculator,  restarting  it  and  reloading  the  program. 


*  This  tape  strip  can  be  obtained  by  contacting  the  Federal  Aviation 
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After  the  primary  register  has  been  loaded  the  secondary  register  is 
loaded  by  depressing  the  following  keys:  "f"  "P-S"*  300  STO  0  200  STO 
4  16  STO  6  1  CHS**  STO  7  1.3  STO  8  "f"  "P-S*. 

Program  Execution 

The  program  is  started  by  depressing  the  "E"  key  for  the  assumed  "E" 
stability  (The  "B",  "C",  and  "D"  keys  will  start  the  program  for  "B", 

"C",  and  "D"  stability  classes  respectively).  The  resulting  three 
numbers  printed  out  on  the  HP-97  or  displayed  on  the  HP-67  after  each 
iteration  is  completed  are;  (1)  time  (in  seconds  )  from  the  program 
start;  (2)  distance  (in  meters)  that  the  point  in  the  emission  tail  has 
moved  down  the  runway;  and  (3)  total  dose  (in  ppm-seconds)  that  the 
receptor  has  received. 

It  is  noted  that  after  20  iterations,  the  dose  value  will  reach  a 
maximum  of  82.16  ppm-sec.  This  value  represents  the  dose  received  at  the 
receptor  from  the  first  emission  release  point  in  the  "tail".  When  the 
dose  converges  on  this  maximum  value  (when  all  concentration  digits 
remained  unchanged  out  to  the  second  decimal  point)  the  R/S  key  is 
depressed  to  stop  the  program.  The  "A"  key  is  then  depressed  to  clear 
registers  and  index  the  starting  point  to  the  second  tail  position.  The 
"E"  key  is  then  depressed  a  second  time  to  start  the  next  computation. 

Again  when  the  dose  value  levels  off  at  35.59  ppm-sec.,  the  "R/S"  key  is 
depressed  to  stop  the  program.  Depressing  the  "A"  key  and  then,  after 
the  display  stops  flashing,  the  "E"  key,  permits  the  last  computation  to 
be  completed  which  results  in  a  dosage  of  40.85  ppm-sec.  for  the  last 
tail  point.  After  reaching  this  last  convergence  value,  the  program  is 
terminated  by  depressing  the  R/S  key.  The  person  making  the  calculation 
can  then  sum  the  three  dose  values  and  divide  them  by  3600  to  produce  a 
one-hour  concentration  of  0.05  ppm. 

CONCLUSIONS 

The  method,  limitations  and  use  of  the  SIMPLEX"A"  model  have  been  described. 
The  program  can  determine  concentrations  from  departing  aircraft  and  has 
the  flexibility  to  easily  accept  parameter  changes.  It  can  treat  either 
single  or  multiple  events  and  permits  air  quality  calculations  to  be 
made  by  persons  without  an  extensive  computer  background.  The  model  can 
assist  in  determining  the  impact  of  aircraft  emissions  on  air  quality  in 
conjunction  with  requirements  for  controlling  engine  emissions  and  can  be 
used  as  a  screening  tool  in  evaluating  the  air  quality  impact  of  proposed 
Federal  actions  at  airports. 


*  The  P-S  command  is  input  by  depressing  the  "CLx"  key  on  the  HP-97  and  the 
"CHI"  key  on  the  HP-67  calculator. 

**  Negative  numbers  are  entered  into  the  Hewlett  Packard  calculators  by 
depressing  the  appropriate  number  key  followed  by  the  "CHS"  key. 


Figure  1 


SOURCE- RECEPTOR  GEOMETRY  DURING  TAKEOFF 


runway 


Case  1 
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Figure  3 

REGISTERS  AND  LABELS 


REGISTERS 
Pr imary 

0  emission  rate 

1  "x"  distance 

2  initial  sigma  "z" 

3  zero  register 

4  sigma  "z"  exponent 

5  sigma  "y"  coefficient 

6  sigma  "y"  exponent 

7  sigma  "z"  coefficient 

8  total  distance  (1+3) 

9  sigma  "z" 

A  sigma  "y" 

B  plume  height 
C  "y"  distance 
D  wind  angle 
E  iteration  time 
I  wind  velocity 

Secondary 

0  fixed  source  receptor  distance  along  runway 

1  plume  rise  factor 

2  dose  summation 

3  sidewind  "y"  factor 

4  fixed  distance  -receptor  to  runway 

5  variable  distance  between  source  and  receptor  in  the  runway 

direction 

6  initial  sigma  "y" 

7  time  at  runway  location 

8  acceleration 

9  hypotenuse  ("h"  in  Figure  2) 

LABELS 

A  program  to  clear  registers  and  make  required  inputs  for 
new  iteration  sequence 

B,C,D,E,  Storage  of  coefficients  and  exponents  for  sigma 
calculations  -  stability  classes  B,C,D,and  E 

4  subroutine  for  Label  A 

5  subroutine  for  Label  A 

6  subroutine  for  Label  A 

7  subroutine  to  switch  registers 

8  main  program  to  move  airplane  along  runway  and  to 

calculate  dosage 

O  ~  C  T  .LmI  Q 
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Figure  4 


INPUTS 


INPUTS  and  OUTPUTS 


Item  Units  Keys 

Primary  Register  (P) 


source  emissions  over 
duration  of  event- 


emission  rate  x  iteration  time 

grams 

Sto  0 

initial  sigma  "z" 

meters 

Stc  2 

plume  height 

meters 

Sto  B 

wind  angle 

degrees 

Sto  D 

iteration  interval 

seconds 

Sto  E 

wind  velocity 

meters  per  second 

Sto  I 

Secondary  Register 


fixed  source  receptor  distance 


along  runway 

meters 

Sto  0 

fixed  distance  from  receptor 

to  runway 

meters 

Sto  4 

initial  sigma  "y” 

meters 

Sto  6 

beginning  time 

seconds 

Sto  7 

acceleration 

meters/ sec/ sec 

Sto  8 

OUTPUTS 

total  elapsed  time  at 

iteration 

seconds 

Prnt  ' 

fixed  source-receptor  distance 

along  runway 

meters 

Prnt  1 

dose  sum 

parts  per  million  sec. 

Prnt  : 

C -lb 


Fi8ure  5 

emission  tail  geometry 


to 


r 
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Figure  6 

SAMPLE  PROBLEM  INPUTS  -  Case  1 


INPUTS 


Item 

Address  Primary  Register 

No./ 

Units 

Keys* 

none 

source  emissions  **over 

duration  of  event 

(emiss.  rate  x  iter,  time) 

(  80  gm/s  x  1  sec.  ) 

80 

grams 

Sto  0 

initial  sigma  "z" 

8 

meters 

Sto  2 

plume  height 

A 

meters 

Sto  B 

wind  angle 

30 

degrees 

Sto  D 

iteration  interval 

1 

second 

Sto  E 

wind  velocity  5 

Address  Secondary  Register 

fixed  source  receptor  distance 

meters  per  second 

Sto  I 

f,  P-! 

along  runway 

+300 

meters 

S  to  0 

fixed  distance  from  receptor 

to  runway 

200 

meters 

Sto  A 

initial  sigma  "y" 

16 

meters 

Sto  6 

beginning  time  *** 

-1 

second 

Sto  7 

acceleration  **** 

1.3 

meters/sec/sec 

Sto  8 

Readdress  Primary  Register  f,  P-S 


*  Applicable  to  both  HP-97  and  HP-67  calculators  except  that  the 

wind  velocity  is  loaded  into  the  HP-67  calculator  by  depressing 
the  black  "h"  key  followed  by  the  black  lettered  "ST  I"  key. 

**  Possible  data  source  -  (7). 

***  For  a  beginning  time  of  zero,  the  negative  value  of  the  iterative 
duration  must  be  input.  This  is  accomplished  by  entering  the 
duration  value  followed  by  the  "CHS"  key. 

****  Possible  data  source  -  Aircraft  Performance  Manuals. 
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SAMPLE  PROBLEM  INPUTS  -  Case  2 


Item 

Address  Primary  Register 

source  emissions  **over 
duration  of  event 
(emiss.  rate  x  iter,  time) 
(  80  gm/s  x  1  sec.  ) 

initial  sigma  "z" 
plume  height 
wind  angle 
iteration  interval 
wind  velocity 

Address  Secondary  Register 


No./  Units  Keys* 

none 


80 

grams 

Sto 

0 

8 

meters 

Sto 

2 

4 

meters 

Sto 

B 

30 

degrees 

Sto 

D 

1 

seconds 

Sto 

E 

5 

meters  per  second 

Sto 

I 

f,  P-S 


fixed  source  receptor  distance 


along  runway 

fixed  distance  from  receptor 

-300 

meters 

Sto 

0 

to  runway 

100 

meters 

Sto 

4 

initial  sigma  "y" 

16 

meters 

Sto 

6 

beginning  time  *** 

-1 

seconds 

Sto 

7 

acceleration  **** 

1.3 

meters/sec/sec 

Sto 

8 

Readdress  Primary  Register 


f,  P-S 


*  Applicable  to  both  HP-97  and  HP-67  calculators  except  that  the 

wind  velocity  is  loaded  into  the  HP-67  calculator  by  depressing 
the  black  "h"  key  followed  by  the  black  lettered  "ST  I"  key. 

**  Possible  data  source  -  (7). 

***  For  a  beginning  time  of  zero,  the  negative  value  of  the  iterative 
duration  must  be  input.  This  is  accomplished  by  entering  the 
duration  value  followed  by  the  "CHS"  key. 

****  Possible  data  source  -  Aircraft  Performance  Manuals. 
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RESULTS  -  Case  1 
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Microcomputer  Graphics  in  Atmospheric  Dispersion  Modeling 


Howard  Segal 

Federal  Aviation  Administration 
WasOington.  DC 


Detailed  atmospheric  dispersion  models  such  as  Airport  Vi 
mty  Air  Pollution  AYAPi  nr  Point  Are.t  l.ine  t PAP '  can 
rijinriusly  mi >del  air  quality  at  man >r  airport.-.  However,  at  the 
>m. tiler  airports  where  fewer  enii-sioti  'ources  are  present  or 
at  the  larger  airports  where  detailed  analyses  may  not  he  re¬ 
quired.  a  well  designed  -v  reening  model  ran  >ave  <  onsiderahle 
time  and  money  Microcomputers  equipped  with  graphics 
tablet'  are  part n  ularly  attractive  for  this  application  -ince 
i.hev  are  inexpensive  ' they  i  an  he  purchased  tor  less  than 
'■'ooio  and  can  in-tuntaneoiisly  accept  -imne  and  receptor 
■ordinates  directly  from  a  base  map  with  little  chance  ol  a 
transcribing  error 

This  paper  describes  the  -pec  lal  leatures  nt  a  multiple 
-oiirce  screening  model  which  ha-  been  programmed  for  a 
mu  roc ompiiter  equipped  with  a  graphics  tablet  and  having 
4yK  of  random  ao.es-  memory  Mk  roc  ••  mi  (niter  i  on  figuration 
i-  -hown  in  Figure  1  The  model  .u  ronym  is  DIMM  iDraphical 
Input  Mic  roci.mputer  Model i  and  it-  -tatus.  features  and 
results  are  presented. 

Model  Status 


Figure  1.  Microcomputer  configuration  tor  ^spe'S'Or.  -ng 


v  eh  ip  emission  and  dispersion  models  for  u>e  hy  field  per¬ 
sonnel 


1 1IMM  i-  an  outgrowth  ot  an  earlier  aircraft  model.  Simplex 
A."  whit  h  was  programmed  tor  a  desk  c  alt  ulator  Simplex 
"A"  was  developed  to  fulfill  the  need  for  a  si  reening  model  lor 
airt  raft  sources.  With  the  advent  of  microcomputers.  Simplex 
A"  wa-  reprogrammed  to  further  simplify  the  assessment 
process.  During  this  reprogramming,  a  method  was  developed 
to  uistantaneouslv  input  source  and  receptor  coordinates  into 
the  microcomputer  directlv  from  a  base  map  with  little  chance 
of  a  transc  rihing  error  This  approach,  which  appeared  soet- 
tee  tive  for  aircraft  sources,  was  then  applied  to  other  airport 
sources  sue  h  as  roadways,  power  plants,  parking  lots,  etc  .  in 
order  to  eliminate  input  errors  and  simplify  model  usage. 
Algorithm'  for  each  of  these  sources  were  then  hatched  to 
■  •  > . . i.a  vck.i.,  a... 


Model  Features 


Point  and  Line  Source  Considerations 


(  oncentrations  from  point  and  line  source-  are  determined 
with  the  classical  point  scan (■  equation  1 


A-  » 


U 

71.1  a  u 


exp 


i  1  . 


The  above  equation,  which  is  used  om  e  tor  e.n  h  point  source, 
is  used  a  number  ot  times  iteratively  lor  line  sources  ('on¬ 
centrations  from  line  sources  are  determined  hv  div  iding  a  line 


Main  terrruna/ 


j  ^  ifg  Departing 

J  y  J  aircraft 


Metro  ra,i 


V—  Aate 


2  Scor:e-re;  ef’Or  aeC^e'r.  ill  Wasl'iccrc'' Nat  2r<H!  ifCOC t 

Area  Source  Considerations 

:r:n  dei  i > j < ■  n  has  not  been  made  oil  whe  ther  to  i  on-idt  r 
-i  •>  1  r>  e-  .I.-  point,'  or  lint-'  Bt  r  ause  point,'  t  an  !*•  pn  s  e»ed 
r  than  line'  ami  lie.  an>r  thtv  can  he  organized  into  am 
eii  area,  a  point  as-umption  wa-  initiallv  u.-ed  Coordi 

■  of  each  [loin!  are  entered  into  the  r  nmputer  v  la  a  pen 
lied  to  ,1  graphic-  Tahir  t  Thi-  operation  i'  quite  tael:  one 

■  r  an  lie  entered  into  the  i  nmputer  and  i  om  duration- 
dated  m  It"  than  three  -er  onds  While  mnt  points  were 
to  'i m (date  the  'ho rt  •  term  parking  area  at  V\  a>hingtnn 
■nai  Airport  1 1 '(  A  >  <  Figure  J  i.  the  ti't'  of  add  it  ional  1  or 
r  i  points  i'  under  tin  r-t  igat  inn  to  help  arrive  at  an  opt  i 

number  of  point-  to  etuplov  in  area  soune  -imiil,i 


■■tnploV'  the  Pillowing  equation. 
■  Vi  ....  [ 


'  -  o  ’•  I  —  I 


The  d.-e-  l rom  e.n  h  1  -  e in: "ion  putt  are  -it  muted  to  give 
the  total  d"'f  a*  a  re.  ept-  r  dm  t..  a  ,  nmplete  take- •((  event 
The  program  -t.  p-  when  the  im  renn  ntal  do-e  increase  he 
.  •  ane-  tn-ignificant  Sim  e  the  air  qualitv  unpai  t  ot  1  he  man  >r 
nreralt  pollutant'  heronte-  m-ignitii  ant  at  onlv  'h. .rt  de¬ 
puties  above  1 1  •.  ground.  '  .  limb  out  and  de-oent  aL'< hm- 
have  been  omitted  !r>  m  the  <  ilMM  program 


Accelerating  Point  Source  Considerations 

IT,*  algorithm  tor  a.  celerat. mt  point  sources  is  described 
del  ail  in  the  Simplex  "A”  t  -er~  <  unde  I  his  algorithm  is 
-ed  upon  the  assumption  that  an  accelerating  point  -mine 
e  .  a  taking  oil  air.  raft  >  relea-t-  it-  emission-  as  a  -ene-  ot 
,-e-  Figure  :  -how-  the  -otin  e  reeeptor  geometrv  whnh 


Receptor  ^ 


In  order  t. .  determine  how  Ions;  it  would  take  n,  enter  data 
in  t I  MM  and  run  it.  a  -tandardi/ed  -r  enario  h-ted  in  the  l’AI, 
l  -er  -  f  i mde  was  u-ed  Thi-  -lenarm  con-Med  ot  a  combi 
nation  ot  I'd  point,  area.  line,  and  a.  .eleratmg  aircralt  -ources 
and  .'<  receptors.  (ilMM  accepted  these  data  in  !’’  ir.'i  and 
call ulated  concentrations  in  l:  _h  Corresponding  times  tor 
an  actual  airport  '  OCA  t  were  |u  min  and  1  ,  h  re-pei  t ivelv. 
Aircraft  sourer—  alone  were  then  run  with  the  (ilMM  mode! 
lor  !  source -reeeptor  scenaru>s  li-teri  in  the  l’AI.  1  -er  -  ( unde 
Two  versions  ul TIMM  were  empluved  the  ! ir-t  with  lower 
•  omputat tonal  ai  ciir.iiv  hut  higher  computer  execut  ion  speed, 
anil  the  second  with  hither  i  omputation.il  accurai  v  hut  lower 
i  omputer  execut  ion  speed  The  higher  -peed  version,  w  hit  h 
increased  exer  ution  -peed  hv  go' . .  was  an  ompltshed  through 
a  power  law  lit  rather  than  a  quadrat ii  lit  to  the  Pasquill 
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Nomenclature 


H  =  effective  height  of  emissions 
Q  =  emission  rate 

Q  i  =  total  emissions  released  during  a  finite  time  period 
u  =  wind  speed 
-V  =  receptor  concentration 

x  =  downwind  distance  in  the  mean  wind  direction 
>  =  crosswind  distance 

n .  =  standard  deviation  of  plume  concentration  in  the 
crosswind  direction 

o  =  standard  deviation  of  plume  concentration  in  the 
vertical  direction 
t 'r  -  receptor  exposure  or  dose 
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The  air  quality  inpact  of  the  proposed  runway  expansion  program  at  Stapleton 
International  Airport  is  determined  in  this  report.  The  method  of  analysis 
is  to  model  the  dispersion  of  pollutants  from  motor  vehicles  and  aircraft 
under  both  1-hour  and  8-hour  worst  case  conditions. 

Results  show  that  aircraft  pollution  concentrations  are  reduced  and  in  some 
cases  completely  disappear  when  the  new  runways  are  added.  This  is  caused 
primarily  by  a  reduction  in  takeoff  delays,  which  are  a  major  objective  of 
the  runway  expansion  program  at  the  airport. 
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SUt^RY 

Most  pollution  -front  aircraft  at  Stapleton  International  Airport  (DEN)  is 
the  result  of  pre-takeoff  delays.  These  delays  result  in  aircraft  queues 
which  increase  the  time  that  aircraft  engines  must  operate  on  the  ground. 
In  just  about  every  case  these  queues,  and  the  pollution  they  create,  are 
reduced  or  completely  eliminated  when  new  runways  are  added. 

This  conclusion  was  based  upon  current  estimates  of  peak  hour  motor 
vehicle  and  aircraft  activity  at  the  airport  and  the  application  of  these 
peak  hour  values  to  all  hours  modeled.  This  approach  provided 
conservative  air  quality  estimates  for  the  two  8-hour  me teorol g i cal  data 
sets  provided  by  the  Colorado  Department  of  Health  (CDH)  since  both  data 
sets  extended  into  the  late  evening  hours  where  there  was  little  aircraft 
or  motor  vehicle  activity. 

A  major  result  expected  from  the  runway  expansion  at  Stapleton 
Internat i on  a  1  Airport  is  reduced  delays  and  therefore  reduced  pollution 
from  aircraft. 


BACKGROUND 

At  the  request  of  Colorado  Department  of  Health  (CDH),  the  Federal 
Aviation  Actni n i strat i on  (FAA)  conducted  an  analysis  of  motor  vehicle 
pollution  at  DEN.  This  study  was  performed  in  conjunction  with  the  runway 
expansion  investigations  at  the  airport.  On  July  1,  the  documentation  of 
this  study  was  transmitted  to  the  CDH  as  report  FAA-EE-86-7  (Reference  1). 

After  reviewing  this  report,  the  CDH  requested  the  following  information: 

1.  An  assessment  of  pollution  from  both  aircraft  and  motor  vehicles 
using  two  sets  of  8-hour  meteorological  data  provided  by  the  CDH. 

2.  An  assessment  of  the  pollution  from  motor  vehicles  and  aircraft 
using  1-hour  ’worst  case*  meteorology.  This  assessment  would 
consist  of  the  addition  of  aircraft  to  the  motor  vehicle  analysis 
of  Reference  1  . 

3.  The  expansion  of  the  motor  vehicle  analysis  of  Reference  1  to 
include  a  wider  variety  of  wind  directions.  The  CDH  recommended 
that  wind  directions  of  180,  200,  225  and  330  degrees  be  modeled 
since  only  westerly  wind  directions  (240  and  270  degrees)  were 
modeled  in  the  original  study. 

4.  Nitrogen  oxides  (N0X>  estimates. 

5.  A  determination  of  the  air  quality  impact  of  motor  vehicles  at 
_ the  1-70/Quebec  Street  interchange.  (Vehicular  flow  rates  to  be 


E-6 


DISCUSSION  -  . 

The  air  quality  impact  of  the  runway  expansion  program  at  DEN  was 
determined  by  calculating  carbon  monoxide  <C0)  concentrations  at  seven 
receptors  placed  in  the  terminal  area.  The  geometrical  location  of  these 
receptors  as  related  to  the  runways  and  roadways  at  the  airport  is  shown 
in  Figures  1  through  5.  The  tool  used  in  assessing  pollution  at  the 
airport  was  the  Graphical  Input  Mi crocompu ter  Model  (GlhM)  (Reference  2). 

Two  "worst  case"  scenarios  were  prepared  in  order  to  calculate 
concentrations  from  aircraft  and  motor  vehicles  operating  at  the  existing 
and  expanded  runway  systems.  Results  from  the  first  scenario  analysis, 
which  included  weather  observations  for  two  specific  8-hour  time  periods, 
are  plotted  in  Figures  6  and  7.  In  all  instances  aircraft  concentrations 
were  reduced  when  the  new  runways  were  added. 

Results  from  the  second  scenario  analysis,  which  employed  estimated  "worst 
case"  1-hour  meteorology,  is  plotted  in  Figures  8  through  21.  The  highest 
combined  aircraft-motor  vehicle  concentrations  are  shown  in  Figure  17. 
These  concentrations  are  significantly  reduced  with  the  introduction  of 
the  expanded  runway  system.  Appendix  B  lists  the  computer  reports  from 
which  Figures  6  through  21  were  prepared. 


MODEL 

GIMM  is  a  complex  source  em i ss i ons/d i sper s i on  model  with  an  emissions 
front  end  which  allows  fast  and  accurate  data  entry  and  "what  if" 
analysis.  The  model,  described  in  detail  in  Reference  2,  is  conceptually 
displayed  in  Figure  22. 

Glhtl  is  compared  to  two  Environmental  Protection  Agency  (EPA)  models: 

Po i n t-Area-L i ne  (PAL)  and  H1UAY  2  in  Figure  23  and  Reference  2.  The 
comparison  shows  Glhtl  results  to  be  very  close  to  those  of  the  EPA 
mode  1 s . 

Before  running  GIMM  it  was  necessary  to:  (1)  establish  source  and 
receptor  locations,  <2)  estimate  vehicular  activity,  and  (3)  select  source 
emission  rates. 

Source  and  receptor  locations  are  shown  in  Figures  1  through  5,  and  the 
rationale  for  developing  data  on  aircraft  and  motor  vehicle  activity  and 
emission  rates  is  described  below. 


DATA  DEVELOPMENT  -  AIRCRAFT 


The  EPA  has  identified  four  operational  modes  for  aircraft  pollution 

aecAceman  f  nnnnneae  •  i  al/  ami  1  r  1  imKnnt  . .**».*!  4.1  * . 


Only  the  queue  and  takeoff  modes  are  included  in  the  model  analysis 
because  climbout  and  approach  contribute  very  little  to  the  pollution 
burden  at  an  airport  (Reference  3).  The  queue  times  selected  for  this 
study  were  15  minutes  for  the  existing  runways  and  3  minutes  for  the 
proposed  runway  configuration.  These  times  are  consistent  with 
capacity/demand  estimates  in  Reference  6.  Appendix  C  describes  the 
adaptation  of  Reference  6  data  to  this  study. 

UJhen  making  a  screening  analysis,  conservative  estimates  should  be  used. 
Peak  vehicular  activity  was  therefore  used  throughout  this  study.  The 
peak  hour  activity  of  aircraft  was  determined  after  reviewing 
documentation  on  actual  aircraft  departures  and  estimates  of  these 
departures  that  were  listed  in  computer  printouts  from  the  Official 
Airline  Guide  (OAG).  Aircraft  activity  at  1700  hours  on  August  19,  1984, 
was  selected  for  this  study.  Eighty-one  commercial,  general  aviation,  and 
air  taxi  aircraft  were  estimated  to  depart  from  DEN  during  that  hour. 

Emission  rates  were  extracted  from  Reference  4. 


DATA  DEVELOPMENT  -  MOTOR  VEHICLES 

Roadway  activity  was  calculated  from  hourly  traffic  counts  obtained  from 
Centennial  Engineering  Company.  Two  traffic  count  data  sets  were  provided 
by  Centennial  Engineering — one  covering  motor  vehicle  activity  on  city 
streets  and  the  other  covering  activity  along  terminal  roadways.  Parking 
lot  activity  was  also  observed  on  August  9,  1986.  Traffic  counts  at  the 
Quebec/Interstate  70  (1-70)  interchange  were  also  provided  verbally  by 
Centennial  Engineering.  From  these  data,  a  roadway  throughput  analysis 
was  prepared  and  vehicular  flow  on  each  roadway  segment  was  determined. 

The  results  are  listed  in  Figure  5. 

Traffic  counts  on  December  20,  1985,  and  August  9,  1986--two  peak  activity 
times  for  motor  veh i c 1 es--were  used  for  this  roadway  analysis.  The 
traffic  at  the  I-70/Quebec  interchange  was  not  included  in  the  modeling 
analysis  since  we  had  not  received  these  data  from  Colorado  State 
personnel  at  the  time  the  model  was  run.  This  information  can  be  easily 
added  to  the  study  when  received.  Because  of  the  great  distance  between 
this  interchange  and  the  terminal,  concentrations  should  change  little 
when  this  additional  data  become  available. 

Peak  hour  activity  was  assumed  for  all  hours  modeled  regardless  of  whether 
the  hour  modeled  was  at  a  peak  ualue  or  not.  This  epproach  provided 
conservative  air  quality  estimates  for  the  8-hour  data  sets  since  both 
8-hour  data  sets  extended  into  the  late  evening  hours  when  there  is 
significantly  less  than  peak  hour  activity. 

Emission  rates  are  calculated  by  a  Mobile  3  submodel  of  611*1.  The 


Item  1  —  An  assessment  of  pollution  from  both  aircraft  and  motor 
vehicles  using  two  sets  of  8  consecutive  hours  of  meteorological  data 
provided  by  the  CDH. 

Most  pollution  from  aircraft  is  the  result  of  pre-takeoff  delays. 

These  delays  result  in  aircraft  queues  which  increase  the  time  that 
aircraft  engines  must  operate  on  the  ground.  In  just  about  every 
case,  these  queues  and  the  pollution  they  create,  are  reduced  or 
completely  eliminated  when  new  runways  are  added.  A  major  result  of 
runway  expansion  at  Stapleton  International  Airport  will  be  reduced 
delays  and  therefore  reduced  pollution  from  aircraft. 

The  air  quality  impact  of  motor  vehicles  alone  is  documented 
in  Reference  1  and  Item  3  below. 

Item  2  —  An  assessment  of  the  pollution  from  motor  vehicles  and  aircraft 
using  1-hour  "worst  case"  meteorology.  This  assessment  would  consist  of 
the  addition  of  aircraft  to  the  motor  vehicle  analysis  of  Reference  1. 

The  conclusion  of  Item  1,  which  was  for  the  8-hour  analysis,  also 
applies  to  the  1-hour  analysis. 

Item  3  —  The  expansion  of  the  motor  vehicle  analysis  of  Reference  1  to 
include  a  wider  variety  of  wind  directions.  The  CDH  suggested  the 
modeling  of  180,  200,  225,  and  330  degree  wind  directions  because  only 
westerly  wind  directions  (2A0  and  270  degrees)  were  modeled  in  the 
original  study. 

After  modeling  the  dispersion  of  pollutants  under  the  four  additional 
meteorological  cases  noted  above,  the  assumption  in  Reference  1  that 
the  highest  concentrations  would  occur  at  the  three  receptors  closest 
to  the  terminal  was  confirmed.  However,  the  wind  angle  at  which  peak 
concentrations  occured  changed.  Revised  peak  concentrations  are  as 
follows:  Receptor  1-30  mg/m3  at  a  wind  angle  of  330  degrees; 
Receptor  2-30  mg/m3  at  a  wind  angle  of  200  degrees;  and  Receptor  3  - 
29  mg/m3  at  a  wind  angle  of  240  degrees.  These  values  were  obtained 
from  the  passenger  vehicle  listings  of  Appendix  A  plus  an  assumed 
average  concentration  for  buses  of  1  mg/m3. 

Item  4  —  NOx  Estimates 

For  all  the  modeling  runs,  NOx  as  well  as  CO  concentrations  were 
printed  out.  Appendix  B  lists  these  data. 

Item  5  —  A  Determination  of  the  Air  Quality  Impact  of  Motor  Vehicles  at 
the  1-70  Quebec  Street  Interchange.  (Vehicular  flow  rates  to  be  provided 
by  the  state.) 

This  portion  was  not  completed  because  traffic  counts  which  were  to  be 
provided  by  Colorado  State  personnel  were  not  received.  However,  it 


TAKEOFF  RUNWAYS  AND  QUEUE  LENr.THS  -  STAPLETON  INTERNATIONAL  ATRPORT 


AKEOFF  RUNWAYS  AND  HUEtlF  LENGTHS  -  STAPLETON  INTERNATIONAL  ATRPORT 
north  departures  -  expanded  runway  system  -  peak  aircraft  activity 


RUNWAY 


TAKEOFF  RUNWAYS  AND  QUEUE  LENCTIIS  -  STAPLETON  INTERNATIONAL  AIRPORT 


RUNWAY 


TAKEOFF  RUNWAYS  AND  QUEUE  LENGTHS  _  STAPLETON  INTER NATT  ON AT,  ATRPORT 
(east  departures  -  expanded  runway  system  -  peak  aircraft  activity) 


TERMINAL  ACCESS  ROADWAYS  -  STAPLETON  INTERNATIONAL  AIRPORT 

(view  A-A  of  Figures  1-4) 
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Peak  hour  activity  assumed  for  all 
hours  of  the  day  modeled 
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APPENDIX  A 


PRINTOUT  OF  KOTOR  VEHICLE  POLLUTION 
FOR  DIFFERENT  KIND  DIRECTIONS  -  APPENDAGE 
TO  REPORT  FAA-EE-86-7  (REF.  I) 
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This  appendix  is  a  printout  of  the  motor 
vehicle  pollution  for  roadways  at  Stapleton 
International  Airport.  It  repeats  the  printout 
of  wind  directions  240  degrees  and  270  degrees 
in  Ref.  1  and  adds  the  printouts  for  wind 
directions  of  180,  200,  225,  and  330  degrees. 
The  date  header  on  each  printout  represents 
the  day  that  the  run  was  made. 
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TYPICAL  PRINTOUT  (THERE  ARE  6  PRINTOUTS  IN  THE  COMPLETE  APPENDIX  WHICH  IS 

INCLUDED  IN  REPORT  EAA-EE-86-11 


P.iSSIK  RATE': 


^ ■ 
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APPENDIX  B 

GIKM  PRINTOUTS  FOR  THE  STAPLETON 
INTERNATIONAL  AIRPORT  SCENARIO 
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Figures  6  through  21  were  prepared  from  the 
printouts  listed  in  this  Appendix.  These  data 
include  the  following  combinations  of  information: 

1.  8-hour  and  1-hour  analysis 

2.  motor  vehicles  and  airplanes 

3.  existing  and  future  runways 

4.  north  and  east  departures 

The  interrelationship  of  these  combinations  is  shown 
in  Figure  B-l . 

To  facilitate  Graphics  Tablet  use,  the  wind  direction 
values  listed  in  the  printouts  had  to  be  referenced 
from  the  top  of  the  page.  This  required  rotating 
the  maps  of  Figures  1-5  from  a  vertical  north 
orientation.  The  user  of  Appendix  B  data  must 
therefore  subtract  90  degrees  from  all  listed  wind 
angles  to  establish  the  true  wind  angle. 

The  graphics  tablet  relates  all  coordinates  to  a  (0,0) 
map  origin.  Since  the  origin  of  the  large  scale  map 
(Figure  1)  is  (0,0),  and  aircraft  coordinates 
are  entered  from  this  man,  these  sources  do  not  have 
to  be  corrected.  However,  the  coordinate  printouts 
for  motor  vehicles  have  to  be  corrected  because  thev 
are  entered  into  the  Graphics  Tablet  from  the 
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COMPUTER  PRINTOUT  SUMMARY 
8  Hour  Analyses 
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Existing  and  Future  Runways 

Existing  Runways 

Future  Runways 

East  Departure 

Run  #5 

Run  #1 

Run  #2 

North  Departure 

| 

Run  #5 

i 

! 

1 

■ 

Run  #3 

Run  #4 

1  Hour  Worst  Cast  Analsvis 


MOTOR  VEHICLES 


AIRPLANES 


Existing  and  Future  Runways 


Existing  Runways _ Future  Runways 


East  Departure 
Run  #10 


Run  #6 


Run  { 1 


North  Departure 


Run  #8 
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TYPICAL  PRTNTOV'  (THERE  ARE  OYER  100  OF  THEM  IN  THE  COMPLETE  APPENDIX) 
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1.  Introduction 


The  aircraft  pollution  burden  is  determined  by  modeling  aircraft  queuing 
emissions.  A  capacity/demand  analysis  is  required  to  estimate  this 
burden.  The  demand  portion  of  this  analysis  is  described  in  section  3.3.1 
of  the  main  report.  The  capacity  analysis  is  made  by  selecting 
north-south  and  east-west  runway  configurations  from  Reference  6  that  were 
determined  to  be  "worst  case"  with  respect  to  air  quality  and  therefore 
appropriate  for  air  quality  analysis. 

2.  Runway  Usage  Selection  for  Air  Quality  Analysis 

Departures  to  the  east  or  to  the  north  were  selected  as  the  approoriate 
runway  geometries  to  use  in  this  air  quality  analysis.  These  two 
geometries,  which  were  part  of  the  seven  geometries  listed  in  Reference  6, 
were  selected  because  they  place  aircraft  emissions  closest  to  the 
pollution  receptors  at  the  terminal  and,  therefore,  would  record  the 
highest  possible  pollution  values. 

For  each  of  the  runway  geometries  listed  above,  two  capacity  values  are 
listed  in  Reference  6;  one  during  flight  under  Visual  Flight  Rules  (VTR) , 
and  the  other  during  flight  under  Instrument  Flight  Rules  (IFF.)  .  The 
values  listed  under  VFR  were  selected  for  this  analysis  because  they 
predominate  under  the  "worst  case"  meteorological  conditions  provided  by 
the  CDPH. 

3.  Delay  calculations 

As  was  mentioned  in  the  main  text  of  this  report  (section  3.3.1), 

81  aircraft  were  estimated  to  depart  from  the  airport  during  the  peak 
hour.  For  the  existing  configuration,  these  aircraft  are  assumed  to 
depart  on  two  parallel  runways  either  to  the  east  or  to  the  north. 
Departures  would  therefore  consist  of  40  airplanes  on  one  runway  and  41  on 
the  other. 

During  the  peak  hour  there  would  be  pressure  to  disperse  these  aircraft  to 
the  third  proposed  runway.  Assuming  an  even  split  over  the  three  runways, 
the  departure  rates  would  be  27  aircraft  per  hour  on  each  runway. 

The  total  VFR  capacity  from  Reference  6  is  150  aircraft  per  hour. 

Assuming  the  takeoff  portion  of  this  capacity  is  88  aircraft  per  hour  or 
9lightly  greater  than  one  half,  the  following  equation  from  Reference  7 
can  be  used  to  calculate  delays  prior  to  takeoff: 

T = 


%  (&0  *r\  i  A,} 

Q(Q~  l) 


q  -  demand  (aire./hr) 

Q  -  capacity  (airc/hr) 
T  »  queue  time  (min.) 
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DELAY  CALCULATIONS  TOR  EXISTING  RUNWAY  SYSTEM 

Assuming  a  Poisson  distribution  of  aircraft  arriving  at  aach  of  the  two 
takeoff  queue  areas: 

T  4-0  (  fe  O  *1  irt . )  .  / 

4-4-  —  4-0)  ~  ^  >n'  (runway  1 ^ 

-r_  4-i  C&>0  O  .  , 

I  ~  ^  Z  4.  (  y  —  I  1  >r'»r'  /  &  I  RC .  (runway  2  } 


For  conservatism  a  15  minute  queue  time  was  selected.  With  a  departure 
rate  of  1  1/2  airplanes  per  minute  the  peak  queue  length  would  be 
10  airplanes.  This  value  is  consistent  with  queue  lengths  reported  by 
tover  personnel  during  peak  hours.  It  is  now  possible  to  UBe  this  takeoff 
capacity  estimate  to  calculate  the  decrease  in  queue  time  with  the 
additional  runways. 

DELAY  CALCULATIONS  FOR  EXISTING  +  PROPOSED  RUNWAYS 

Assume  that  the  proposed  east-vest  or  north-south  runways  are  in  place  and 
that  during  peak  hours  scheduled  departures  will  be  evenly  directed  to 
these  three  runways  (27-27-27).  Assume  that  each  runway  has  a  takeoff 
capacity  of  44  (1/2  of  88)  departures  per  hour. 


Under  these  conditions: 

“T  2  1  (  ) 

'  ~  4_4-  C^-ZT)  ~  Z  l(> 


(all  3  runways  ) 


To  be  conservative  3  minutes  was  selected 
Therefore,  aircraft  queues  at  each  runway  will  be: 

PRESENT  CONFIGURATION  -  15  minutes 

PROPOSED  CONFIGURATION  -  3  minutes 
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